Self-sustained Clusters and Ergodicity Breaking in Spin Models 
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Tlie emergence of self-sustained clusters and their role in ergodicity breaking is investigated in 
fully connected Ising and Sherrington-Kirkpatick (SK) models. The analysis reveals a clustering 
behavior at various parameter regimes, as well as yet unobserved phenomena such as the absence of 
non-trivial clusters in the Ising ferromagnetic and paramagnetic regimes, the formation of restricted 
spin clusters in SK spin glass and a first order phase transition in cluster sizes in the SK ferromagnet. 
The method could be adapted to investigate other spin models. 
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Spin glasses are magnetic materials characterized by 
extremely slow magnetization relaxation in the absence 
of external field [l|, 0] . Several models have been devel- 
oped to explain their behavior 0, 3 which in turn have 
revealed a rich physical picture of a rugged free energy 
landscape 0, Q ■ Remarkably, the physics of spin glasses 
has a non-trivial connection to interdisciplinary applica- 
tions including image processing, error correcting codes, 
neural networks and combinatorial optimization [a, 0|. 
Its connection to structural glasses and supercooled liq- 
uids have also been explored to explain the physics below 
the glassy temperature p, M, ■ 

Among the various spin glass models, the Sherrington- 
Kirkpatrick (SK) model 0] is arguably the most studied. 
One of the most intriguing features of large scale disor- 
dered systems in general and the SK model in particular, 
is the breaking of ergodicity in some parameter regimes 
(e.g., temperature, strength of interactions or topology), 
particularly in the spin-glass phase where it manifests it- 
self through a complex symmetry structure of order pa- 
rameters that describe macroscopically the correspond- 
ing solution space. Although macroscopic properties of 
the SK model are relatively clear its microscopic features 
are less understood [lOJ, in particular the existence of 
stable domains that are independent of the remainder of 
the system; these are important for gaining insight into 
the mechanism that gives rise to ergodicity breaking and 
the physical picture of spin glasses. 

In this Letter, we examine analytically the existence of 
self-sustained spin clusters in fully connected Ising and 
SK models. We remark that a similar behavior, termed as 
backbones or frozen variables in sparse systems [lll.ll2|. is 
induced by the topological disorder and is therefore some- 
what different from the self-sustained clusters studied 
here; nevertheless, sparse topologies could be analyzed 
by extending the method presented here. We study the 
existence and nature of self-sustained clusters in various 
phases, the dependence of their sizes on system param- 
eters and the existence of phase transitions with respect 
to cluster sizes. The Ising model will be analyzed first, 
followed by a more involved analysis of the SK model. 

Models - The SK model comprises TV spin variables, 
any two of which i and j interact via a ferromagnetic 
{Jij > 0) or anti-ferromagnetic ( J^ < 0) symmetric cou- 
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FIG. 1: Phase diagram of the SK model as a function of 
coupling mean Jo and temperature T. Right hand side of 
the dashed line corresponds to the region where a first order 
phase transition in cluster sizes is observed, i.e. self-sustained 
clusters of certain size are absent. 



pling {Jij — Jji). Coupling variables are randomly drawn 
from a Gaussian distribution of mean Jq/N and vari- 
ance J"^ /N] the corresponding Hamiltonian is given by 
"HsK = ~S(i7) JijSiSj, which sums over all un-ordered 
spin pairs {i ^ j). The infinite-range Ising model is a spe- 
cial case of the SK model with J = or Jo — > oo, and the 
corresponding Hamiltonian is Hising= ^Jo^a-]) SiSj/N. 
To compute physical quantities of interest, one uses 
the identity IwZ = lim„_>.o(-Z" — 1)/'^ to carry out an 
averages over quenched variables, replacing the average 
of In Z by that of the replicated partition function Z"^ [5| . 
As A^ ^^ oo, solutions space is described by the magneti- 
zation and inter-replica spin correlation order parameters 



^2^SiQ., qap- ^l^Sia 



Si/S, 



(1) 



where a, /3 = 1, . . . , n are replica indices. An ultramet- 
ric structure of the order parameter symmetry is then 
used to facilitate the calculation, the simplest of which 
is the replica- symmetric (RS) ansatz, where one substi- 
tutes rria = m for all a and qa/j = q for all a ^ (3. The 
various phases observed in the model are expressed by 
the values of m and q, for instance the paramagnetic 
{m = q = 0), ferromagnetic (m ^ 0,q> 0) and spin glass 
phases {m = 0,q>0) as shown in Fig. [ij 

Self- sustained clusters- Denote a set C of spin variables; 
for each spin i G C we define in-cluster and out-cluster 



magnetic fields Ui ~ J2jeC '^ij^j ^^'^ '^i ~ '^j^c '^v^j 
induced by spins in and out of C, respectively. The total 
magnetic field experienced by spin i is hi ~ ui + Vi. The 
set C is self- sustained if 



> \V^ 



VieC. 



(2) 



In other words, the magnetic field experienced by each 
individual spin i in C is dominated by the contributions of 
peer spin variables in C. We remark that our framework 
can accommodate other cluster definitions. 

To obtain the distribution of clusters, we denote r2(r) 
to be the number of self-sustained clusters of normalized 
size r = \C\/N. Since In ri(r) is an extensive quantity, we 
define the entropy of clusters to be S{r) = [lnO(r)]/A^. 
For instance, one can easily compute S{r) of the Ising 
model at zero temperature T = where all spins are 
aligned. Since the couplings are uniform, Eq. ([2|) is 
satisfied for a set C if r > 0.5. Indeed, any grouping 
with at least half of the spins is self-sustained, which 
implies n{r) = C^^ ^ N\l[{rN)\{N - riV)!] and S{r) == 
— r In r — (1 — r) ln(l — r) for r > 0.5; and ^{r) =0 and 
S{r) — — cxD otherwise. We note that using this defini- 
tion, self-sustained clusters which arc subsets of larger 
self-sustained clusters are also counted. 

We further define a variable ai =- 1 , — 1 to identify cases 
when spin i is included in or excluded from the cluster, 
respectively. Thus, the cluster size 'r = ^^(\+a.i)l2. One 
can then define an indicator function 



'({a,},{sJ,{J.,})-n 



1-0-, 1- 



^e 



'v^-v?) 



(3) 



where the step function Q{x) = 0, 1 for the cases .t < 
and X > 0, respectively. It turns out that the value of 0(0) 
is crucial in the paramagnetic phase as will be discussed 
later. Thus, ui = 1 if the cluster defined by the set {ai — 1} 
is self-sustained, and w = otherwise. 

Ising Model - To derive S{r) for the fully connected 
Ising model at any temperature T we uniformly sample 
spin configurations of given magnetization fn = '^^ Si/N , 
as it uniquely defines the model's macroscopic properties. 
It is sufficient to introduce an operator partition function 
which measures the entropy S{r) of clusters given m: 



Rising (7,™) 



(4) 



N 



.. Ei(i^i) 



where the dependence of w on {Jij} is omitted as they 
are all identical (Jo)- The parameter 7 plays the role of 
pesudo-temperature conjugate to the cluster size X]i(l + 
cri)/2; by computing Z^ one obtains the entropy £'(7) and 
cluster size r(7) as a function of 7 leading to S{r). 

Details of the calculation arc found in the Supplemen- 
tary Information (SI); here wc briefiy describe the solu- 
tion. In the limit of A^ — ^ oo, Z\sing is given by 



-Zising(7, m) = yl(TO)[l -I- Q{msam)e^]^ 



(5) 



where the variable rusa 

given sclf-consistently by the equation 

f Q{mscrm)e^ - 1 

nisa ~ rn -—, 

\Q{m.sam)e^ + 1 

The prcfactor Aim) in Eq. ([5]) is given by 

A{m) = e-'^^'^°™'[2cosh(/3Jom)]^ 



X^i SiUi/N and its value is 

(6) 

(7) 



the entropic contribution of spin configurations {si\, i.e., 
A{m) = Tr{si} ^{^iSi/N — m). Indeed, the partition 
function of the Ising model is e~^-^'='"K^(?Ti), with average 
energy £;isi„g = -NJom'^/2. 

Using Eq. ([5]), one can drive the cluster size r{j) by 



nsineh,m) = 



1 dhlZisinn 



N 



d"f 



Q{ras„m)e~' 
Q{msam)e'^ -f 1 



(8) 



To compute the entropy 5(7), one subtracts the entropic 
contribution In^ from In^^isjng and apply the Legendre 
transformation to obtain 



Slsingil.m) 



1 

N 



7 — — ^ +ln Zising - In A 



7 &{msam)e"' 



(9) 



hi[l + Q{m 3^^11)6^ 



e(m^<,m)eT + 1 

To obtain S{r), we assume m>Q and solve Eq. (jH]) to 
obtain irisa = —rn, for 7 < 0, mga = rn ( 1-,^\ ) for 7 > 



e-i+l 



and no solution in the range —mKrag^ <0; as a result 







S{r) 



r = 

-00 < r < 0.5 , (10) 

— r Inr — (1 — r) ln(l — r) r > 0.5 



shown by the black line in Fig. [2] This result is valid for 
the entire ferromagnetic phase {m ^ 0) and is consistent 
with S{r) at r = obtained by simple counting. It 
implies that in the ferromagnetic phase, regardless of T 
and ?n, clusters that include at least half of the spins are 
self-sustained and the magnetization is uniform over any 
subset of spins even for small m values. Alternatively, one 
calculates the in-cluster and out-cluster magnetization^ 



{S.)c 



2r 



2(1 -r) 



, (11) 



respectively, to show [using Eqs. ^ and ([H])] that self- 
sustained clusters have the same magnetization as the 
out-cluster spins, {si)a-i=i = {si)ai=~i—fn- 

For the paramagnetic phase, S{r) is ambiguous since 
TO = and |u,;| = \vi\ = in Eq. ([2]); it thus depends 
on the definition of 9(0) in Eqs. © and ©. Only 
trivial self-sustained clusters are observed: the choice 
0(0) = 1 results in S{r) = — r Inr — (1 — r) ln(l — r) for 
all cluster sizes < ^ < 1, implying that any subset 
of spins is considered self-sustained, while for 0(0) = 0, 
S(fi) = and S{r) = —00, Vr 7^ 0, implying that no 



self-sustained clusters exist. We note that from Eq. (|TT|) . 
{si)ai=i = (si)o-i=-i = regardless of the value of 6(0), 
which implies that magnetized domains are always ab- 
sent from the paramagnetic phase. 

SK Model - Similarly, in the SK model we uniformly 
draw system configurations from a distribution defined 
by the order parameters {ma} and {gap}, and introduce 
an operator partition function which measures S{r) given 
{771a} and {qa/s}- Unlike the Ising model with a single or- 
der parameter m, the order parameters in the SK model 
are labeled by replica indices, we thus define a replicated 
operator partition function 



EsKil.imaj^iqapj.n) 
= TV Tr Tr e'^^-.- 



(12) 
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We further define the corresponding un-replicated par- 
tition function with respect to spin configurations as 
2sKh,P{ma),P{qai3)], such that P{ma) and Piqap) are 
the distributions of nia and qap in the limit n — ?► 0. The 
logarithm In 2^gK is given by 



In ZsK [j,P{ma),P{qal3 )] 



lim 



■^■SK 



(13) 



To find the exact form of P{ma) and P{qap) in the spin 
glass phase requires the full replica sysmetric breaking 
(fuU-RSB) ansatz, which is in principle feasible but very 
difficult. We will thus compute In^^sK under the replica 
sysmetric (RS) ansatz, where P{ma) ^ S{ma—'in) and 
Piqafs) = S{qaf3-q){l - 6a,b) + S{qafi ~ l)Sa,b such that 
ln2^SK only depends on the variables 7, m and q. 

Even with the RS ansatz, the calculation of In^^sK is 
rather involved. We will thus describe the main rationale 
and results and refer readers to the SI for details. To 
obtain S{r), we compute 7'sk(7,"t-,(z) and S'sk(77'ti,(z) 
by similar equations to Eqs. ^ and (|9]) with InZising 
replaced by In Zgx and In A replaced by the spin entropic 
contribution in the SK model 



lTiB{m,q) = lim — 

n~>0 n 



(14) 






N 



YlS{^l^^-q]-l 



a^l3 



N 



Figure [5] shows S{r) as a function of Jq at T = 0.5. 
Remarkably, in the spin glass phase (e.g., Jo = 0.5) clus- 
ter entropies exhibit a similar general shape to those ob- 
tained by counting in a uniform spin configuration but 
with degrees of freedom reduced (almost exactly) by half 



S{r)', 



-7' Inr — (1 — r) ln(l — 7') 



InC 



N/2 
Nr/2 



N 



(15) 



as shown in Fig.|31Ja). We observe that this picture holds 
in the spin glass phase regardless of the values of T and 




FIG. 2: (Color online) The entropy S{r) of self-sustained clus- 
ters of size r at T = 0.5 and various coupling mean values Jo 
in the SK model. Curves with open symbols have been ob- 
tained with a fine resolution of 7 at intervals of 0.02. The 
black line Jo — >■ 00 corresponds to S{r) in Ising ferromagnet. 




FIG. 3: (Color online) (a) Entropy S{r) of self-sustained clus- 
ters in the spin glasse case with Jo = 0.5,r = 0.5, compared 
to the limit of [InC^^ /jJ/A'' as A'' — >■ cxj. Inset: the maximum 
of S{r) obtained by exhaustive search in the ground state of 
small SK systems with Jo — 0.01, compared to the value of 
S{0.5) = llnC^jl]/N = 0.346 (red dashed hne); both mean 
values and standard deviation are shown, (b) The difference 
qscrsa- — qrn^ as a function of r for various Jo values. 



Jq. To test the validity of this result we studied numeri- 
cally S{r) in small SK systems by exhaustive search. The 
inset of Fig. ^a.) shows that max^ [5(r)] is approaching 

[\nCj^,^]/N as N increases in agreement with the theo- 
retical predictions. 

This profile of self-sustained clusters is consistent with 
our understanding of the spin-glass phase: firstly, it 
shows a gap between the trivial cluster that encompasses 
the entire system (r = 1) and the exponential number 
of smaller self-sustained clusters which presumably cor- 
respond to suboptimal solutions; secondly, it shows that 
smaller size self-sustained clusters are determined by ap- 
proximately half of their constituent spins while the other 
half are fixed by inherent system correlations. 

To further understand the relation between self- 
sustained clusters and ergodicity breaking, we examine 
the difference d ~ qsasa — qTi^ where 



qsasa — [{siaO'iaSij3Crii3}i,a,l3], nia = [{(Jia}i 



(16) 



with [. . . ] corresponding to the average over coupling dis- 
orders. One expects d = when the spin-configuration 
overlap between two replica is uncorrelated with cluster 




FIG. 4: (Color online) (a) The self-sustained cluster size rsK 
as a function of 7 at temperature T = 0.5 with Jo from 1.2 
to 1.8. (b) In-cluster and out-cluster magnetizations as a 
function of the cluster size r. 



affiliations; on the other hand, d>{) when correlated spin- 
configuration in two replica tend to have correlated clus- 
ter associations. Figurc|31Jb) shows d as a function of r for 
various Jq values. In the spin glass phase (Jo = 0.5) d>Q 
for all r, suggesting: (i) the presence of self- sustained and 
frozen spin clusters of all size; (ii) not all spin subsets 
constitute self-sustained clusters, or otherwise d would 
have vanished. These results suggest that an extensive 
number of spin flips are required to destabilize or modify 
self-sustained clusters, which points to the existence of 
high energy barriers that lead to meta-stable configura- 
tions. The same phenomenon is identified as backbone 
or rigidity in sparse spin systems studied elsewhere. A 
similar picture emerges in the ferromagnetic phase with 
small Jo > 1, except that (i«0 for small r values. This 
result and the cluster magnetization (si)cri=i — J'O as r^-O 
(see Fig. IH^b)) suggest that small self-sustained clusters 
are not frozen and thus can be easily flipped to merge 
into larger clusters. All these indicate that ergodicity 
breaking is most prominent in the spin glass phase. 

We continue to examine S{r) by increasing Jo, exiting 
the spin glass to the ferromagnetic phase, where one ex- 
pects a different profile of S{r) than that of Eq. (|15p : this 
difference is particularly emphasized when one considers 
the limit of Jo -^ 00, which corresponds to Eq. ([TO)). 
The cluster entropy, shown for increasing Jo values in 
Fig. m exhibits the onset of discontinuity in cluster sizes 
at Jo = 1.8, implying the absence of clusters in a range 
of sizes. The range where discontinuity occurs increases 
with Jo until S{r) reduces to Eq. pH)) when Jq ^ 00 as 
shown in the SI. To examine this behavior we plot the 



expected cluster size r as a function of 7 in Fig. lU^a) 
(higher 7 selects clusters of a larger size). An abrupt 
jump in cluster size appears when Jq > 1.6, resembling 
a first order transition, which implies the emergence of 
large and small clusters and the absence of clusters of 
sizes in between. The phase boundary identifying the 
onset of this first order phase transition was added to 
the SK phase diagram in Fig. [1] denoted by x symbols. 
This phase line marks the emergence of an extensive fer- 
romagnetic domain which grows in size as Jq increases 
and becomes the trivial cluster in the limit Jq— ^00. 

One should note that r(7) is not identically zero be- 
fore the transition point, implying the presence of small 
clusters in the ferromagnetic phase, which presumably 
correspond to small spin domains of arbitrary align- 
ment. Figure IH^b) shows that the in-cluster magneti- 
zation {si)cri=i > ni for the entire range of r except 
when r > 0. This result and the out-cluster magnetiza- 
tion {si)ai=-i <m suggest the presence of local domains 
of weaker magnetic alignment. We remark that similar 
magnetization domains do not appear in the Ising ferro- 
magnet, suggesting that coupling disorder is crucial for 
the formation of such domains. 

Summary - We showed that self-sustained clusters re- 
late to the formation of meta-stable configurations sep- 
arated by an extensive number of variables, one of the 
main features exhibited by disordered systems in the 
spin-glass phase, which leads to ergodicity breaking. 
Such domains have been termed backbone variables else- 
where. We reveal the existence of such clusters in the 
spin-glass and ferromagnetic phases of the SK model and 
the absence of non-trivial clusters in the Ising model. 
Other observations include a first order phase transition 
in the size of self-sustained clusters and the presence of 
domains of stronger magnetic alignment in the SK fer- 
romagnetic regime. The role of self-sustained clusters 
in different spin models, both sparsely and densely con- 
nected, is yet to be investigated analytically for gain- 
ing insights into the corresponding physical behavior; the 
new framework and understanding may also play an im- 
portant role in interdisciplinary applications, particularly 
the development of optimization algorithms. 

This work is supported by the EU project STAMINA 
(FP7-265496) and Royal Society Grant IE110151. 
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O : 1 Derivation of Zi^mg 

!>; '• 

•^ ! We would like to study the properties of self-sustained clusters without affecting the Ising 

system itself. To do that we exploit the fact that the model has been solved previously and 
is determined macroscopically by the order parameter m. We therefore uniformly draw 
configurations s that are consistent with the macroscopic description of the model, the 
value of the parameter m. 

We start from the operator partition function Rising of the Ising model 



^ising(7,M) = Tr Tr ^({a,}, {s,})<5 f ^ - m) e^^^^^, (SI) 

where we denote the magnetization of the Ising model as capital letter M in this supple- 
mentary information, instead of m as in the main paper, to avoid confusion in subsequent 
derivations. The indicator function w is given by 



w iWi}, {si}, Uij}) = n 



l-O-j l + C^^o / 2 2\ 



(S2) 



where the step function 9(x) = 0, 1 for the cases x < and x > respectively, and assume 
6(0) to be either 1 or 0, as discussed in the paper; the variable cxj = —1, 1 corresponds to 
the case when spin i is included in or excluded from the cluster, respectively. The in-cluster 
magnetic field Ui and out-cluster magnetic field Vi are defined by 






iec 



Vi 






(S3) 
(S4) 



m 



Using these definitions, the argument u} — vf in the step function of Eq. (SEJ can be 
simplified as 



2 2 

Ui-Vi 



J^ JijSjaj j f ^ JikSk j 



(S5) 



In order to trace over {sj} and {ai\ in Eq. (31]), one has to factorize the terms in Rising 
over i. To achieve the goal, we: (i) use the integral representation of delta function to 
represent 5(^1,^ Si/N — M); (ii) denote hi = J2j Jij^j ^^^d Vi = Tlj ^ij^j^j- These lead to 



2i 



Isine 



Xr Tr / ^^iMM-iM^f± 



{'^i}W^}J 277 

"■'^jd-'^i / "■Vid-Vi ihjhi-ihi y:, JiiSi+iriifji-ifi, y:, JjiSja, 



X 



X 



n 

i 

n 



277 J 27r 
1 — cTj 1+a 



+ 



e {hir]i 



Si(14<T^) 



gT 2 



(S6) 



We then substitute Jij = Jq/N in the Ising model case and introduce the following mean- 
field parameters 



(S7) 

(S8) 

(S9) 

(SIO) 



rus = 


1 \- 

i 


rusa -- 


i 


^h = 


i 


rrifi = 





One should note that nis is indeed the given magnetization M in Eq. (Sd]), so we expect 
m^ = M to be a consequence of the subsequent derivation. The factor —i included in 



the definitions of rrij^ and rrifj has been introduced to facilitate the calculation later. Using 
again the integral representation of the delta functions for these mean-field parameters, 
Rising becomes 

f dM f dnisdrhs f dmsadrhsfj f druf^drhf^ f dmfjdrhfj 

-^IsiiiK — / ~Z 



^'""^ ' 2tt I 2tt I 2tt I 2tt I 2tt 



X exp \iMM+irnsms+irnsa'tnsa + i'>T^t''n;+i7ni^mi^ — iM^ 



n E E 

1 — (7j 1 + cr, 



dhjdhj I dVi^4i_ ih^hi+iriifii 
2-K / 27r 



X 



+ 



e {hiTji 



■jj[-imsSi-iihsaSiai-im-^(~ihi)-irhfi(-if)i)\+'y^f- 



(Sll) 



With the change of variables M — )■ iNM, rhs — )■ iNrhs, rhsa — > iNrhs^, rhf^ — ;■ iNrh^ and 
rhfi^iNrhfi, one can show that Zismg is given by 



■^Ising '^ 



dM C dmAm^ C dm^„drh 



^sa^"^S(j 



dm-^drhf^ f dmfjdTJifj ^^ 



27r ./ 27r ./ 27r ./ 2n J 2n 

such that 

Rising = —MM — nisihs — ms^ihsa — ^^iJ^h " f^n'ft^fi + Mnis + jQiTisTrij^ + Joms^rrif^ 



(S12) 



+ '°g EE 



dh dh f drj dfj 



'1-cr 1 + 0- 
X 



^ihh+ir)fi 



271 



Q{hT]) 



^msS+ThsaSa+mf^(—ih)+mfj{~ifi)+'y 



i+j- 



(S13) 



where the site index i is omitted as all the terms are factorized. One can then integrate 
h and r) which become the delta functions 6{h — rhfj and 6{t] — rhfi); implementing these 
delta functions by integrating h and r], \E'ising becomes 



Rising = —MM — nisrhs — rris^rhsa — '"^fj^h ~ ^^fj + Mrris + Joirismf^ + J^rris^mf^ 

(S14) 



'og-iEE 

^s=±lo-=±l I- 



1-cr 1 + 0-^,^ „ , 



^msS+msa-S(T+7ii2L 



We can then sum over of s and a, such that ^E'lsmg becomes 

Rising = - MM - nisrhs - nis^rhsa - ^h^u ~ f^v^v + Mnis + Joirismf^ + Joms^rrif, 

+ log [2 cosh(ms — rhsa) + 2 Qifhf^rhf,) cosh(ms + msfj)e'] . (S15) 

To compute 2^ising, one can then make use of Eq. ( 91^ to evaluate Zi^i^g by the method 
of steepest descent, such that the integral is given by e^*^='"e when \E'ising attains its max- 
imum value. We thus differentiate \E'ising in Eq. (S fTSj) with respect to M, the m variables 



•s 


= M + JoTUf^ 


'SO- 


= Jorrifi 


h ^ 


= JoTTls 


V ~ 


= Jorrisa 



and the m variables (but not M which is a given constant). The differentiation of Eq. ( 9T5|) 
with respect to M leads to 

ms = M (S16) 

as expected and suggested by Eq. (311) • The differentiation of Eq. (9T5!) with respect to 
the m variables gives 

(S17) 
(S18) 
(S19) 
(S20) 

The differentiation with respect to the m variables results in 

_ sinh(m^ - rhsa) + Q{rhf/hfi) sinh(m^ + ■msa)e' 

^^s — 1 — TT ~ \ r^ / ^ ~ — \ 1 — T^ ~ \ : V / 

cosh(ms — mso-) + ^yrajnifi) cosh(m5 + msa)^^ 
_ - sinh(mg - m^^) + Q^rhf^rhfj) sinh(mg + msa)e'^ 
cosh(ms — msa) + 'Syrrijjnf,) cosh(ms + nisaje^ 

mf,5{mfjhi^) sinh(ms + msa)e^ 

cosh^m^ — nisa) + 'a[mymfi) cosh^m^ + nisaje'^ 

^_ ^ m^^(m^m^) sinh(m, + m,^)e^ 

^ cosh(m5 — rhsa) + Qifnfrhfi) cosh(ms + ^so-)e'^ ' 

where 6{x) is the delta function. One can summarize all the above relations into 4 equations 
with 4 unknowns, namely M,mscr,fnf^ and rriff, and the following equations 

_ sinh(M + Jpm^ - Jpm^) + 9(Mm^^) sinh(M + Jpm^ + Jom^)e^ 

cosh(M + JoTUf^ — Jomfi) + Q{Mmsa) cosh(M + Jotrif^ + Jomfi)e'^ 

_ - sinh(M + Jpm^ - Jprrifj) + Q{Mmsa) sinh(M + J^m-,^ + Jomfj)e'^ 

cosh(M + Jpwi^ — JoTJirj) + Q{Mmsa) cosh(M + Jpm^ + Jomfj)e"' 

msaS{Mmsa) sinh(M + Jpm/^ + Jpm^)e^ 

m^ = — , ('5^' j 

cosh(M + Jorrif^ — Jpmjj) + Q{Mms„) cosh(M + Jpm^ + Jomfj)e'^ 

_ MSjMms^) sinh(M + Jpnif^ + Jpm^)e^ 

cosh(M + Jom-^ — Jomfj) + Q{Mmsa) cosh(M + Jprn^ + Jpm^)e'^ 

where the original argument of the step function is J^Mmga and we have omitted the factor 
Jp since it is always positive and does not influence the value of the step function. Although 
it seems difficult to solve equations (!25l) - (!28l) . we will show later that m^ = nif, = is a 
self-consistent solution. We thus put m^ = m^j = into Eq. ( 925|l which leads to 

M = tanhM. (S29) 



Since M is the given magnetization of the Ising model, the above equation is satisfied by 

M = /3JoM (S30) 

such that M = tanh(/3JoM) as in the original Ising model, and the physical inverse tem- 
perature /3 appears naturally even if one assumes no knowledge of the temperature T in 
Zising in Eq. (SH]). In this case, from Eq. (^26ll we have 



as stated in the main paper. Finally, we use Eqs. ( 930|) and ( 93T]) to show that the ansatz 
rrij^ = rrifi = is consistent. We note that when M 7^ 0, i.e. in the ferromagnetic phase, 
^^so- 7^ for all 7 except the ambiguity at 7 = — ln[G(0)] which turns out to be the 
singular point at cluster size r = 0.5 in the final solution, thus m^ = m^ oc 6{Mmsu) = 
in the ferromagnetic phase. For the paramagnetic phase, the magnetization M = and 
m^ = rrifi oc ■ 6{0) which do not have a well defined value. Nevertheless, putting m^ = 
nifi = results in Eqs. (S l29|) and (S l30|) as obtained in the original Ising model and is thus 
a consistent solution in the paramagnetic phase. 

We finally substitute Eqs. (9291) -(£131]) and m^ = m^j = into Eq. (SlTSl) which results 
in 

Rising = -(3JoM^ + log{2cosh(/3JoM)[l + e(Mm,,)e^]} . (S32) 

In the limit A^ — ?■ 00, Rising = e^"^^='"s and is given by 

Zi,i„g(7,M) = A{M)[1 + e{Mms.)eY, (S33) 

such that 

A{M) = e-^^-^o*^'[2 cosh(/3JoM)]^ (S34) 

as stated in the main paper. By a similar calculation, one can show that A{M) = 
T^isi} S{J2i ^i/^ ~ M)y i-6- the entropic contribution of the spin variables {si}. 

2 Derivation of Z^k 

To derive 2^sk, we start from the replicated operator partition function 

SsK(7,{Mj,{Q„4,n) = Tr Tr Tr e^E.,. ^-'r'"' J] w({a,J,{s,J,{J,,}) 



i[6 (^^-^^ n^ (^^-g"'^) p(j) 



where we will use the capital letters Ma and Qa^ instead of nia and q'^^ to represent 
the SK model order parameters to avoid confusion in subsequent derivations. Also here, 
we uniformly draw system configurations from those which are consistent with the order 
parameters that fully describe the model macroscopically. Following the expression of 
Eq. (36]) in the case of the Ising model, we: (i) use the integral representation of delta 
function to represent S(^^Sia/N — Mq) and 6 {^^SiaSijs/N — Qa/s); (ii) denote hia = 
Yl,j JijSja and r]ia = Y.j JijSjaCTja- Thcse lead to 
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We then average the coupling disorder by integrating Jij over the distribution 



(S35) 
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for each (ij), to obtain 



(ij) 



I I ^ / rl J ■ • n( J ■ •\ ^ "^^ij Z-^a \ ^ic(^ja.\'^ja.^ia.\'nioL^ja(^ja\'^j(y.^ia.^ia I 



n 



J2 



oc I I exp 



2N 



7 ^ [hiaSjahipSjfi + hiaSjahji3Sij3 + hiaSja'i]il3Sjl30'jl3 + hiaSjafjjpSipa., 



ifi 



{ij) V a,l3 

+ hjaSiahipSjp + hjaSiahjpSip + hjaSiaflipSjpajp + hjaSiafjjpSipaip 
+ fliaSja<^jahii3Sji3 + fliaSja<^jahjf3Sii3 + fjiaS ja<y jofjipS jl3<^ jji + fjiaSja'^jaVjf^SipO'ilS 
''Ija^ia^ia'^il3^jl3 ' 'Ija^ia^ia'''jl3^ifi ' 'IJ0i'^ia^ia'lil3^jl3^jl3 ' 'Ija'^ia^iaVjlS'^ilS^if} 

*"77 / ^ I '''ia^ja ' '''ja^ict ' Via^ja^ja ' Vja^ia^ia j f 
a ) 



iVJ^ 



exp 



a,f3 



NJ< 



^VE[2<K + 2<.K] 



(S36) 



where we have neglected terms of 0{N) in the last line and keep only terms of 0{N'^). We 
then define mean-field parameters, as in Eqs. (S[7|)-(S fT0|) in the case of the Ising model, to 
be 

< = ^E^-' (S37) 

i 
i 

K = -]^E^-' (S39) 



^" = -^E^- (S40) 



" iV 



We also define a set of mean-field parameters to account for the correlation between replica. 



(is 



a0 _ 

N 



^^SiaSip, (S41) 



<r = ^E^-^^/^' (S42) 



^^t = ]7E^-^^/^' (S43) 

i 

9.1 = ^E^-^^/^^^/^' (S44) 
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^^^^ =~iv5Z^^°^^/3' (^^^) 
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a/3 ^ 

^rjsa 



-TT^^Vi^f^ip, (S49) 



5??S(T JY 



rr 2J ■SjaSi/3CriQ(Ti/3 . (S50) 



We note that m" and gf/ are indeed equivalent to Ma and Qa^, respectively, and hence 
we expect to obtain m° = Mq, and g°/ = Qq,^ in the derivation. The above mean-field 
parameters are introduced in Ssk by the integral representation of delta functions, which 



lead to 
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(S51) 



We proceed by: (i) making the change of variables Mq, -^iNMa, Qa/B ^iNQ^p, and 
similar changes of variables for all the other variables of g"^ and rha^ (ii) assuming replica 
symmetry (RS) such that for all a, M^ = M, M^ = M, and similarly for other m^ and rha 
variables, i.e. rria = m and rha = itl; (iii) for a ^ P, Qa/B = Q and Qa^ = Q, and similarly 
for other g"'^ and g"^, i.e. g"'^ = g and g"^ = g; (iv) for a = (3, Qa/B = 1, Qa/3 = C, and for 
other variables of g"^ and g"'^ we assume g"'^ = c and q'^^ = c. In this case, one can show 
that 



.=,SK oc e 



Nn'i'sK 



(S52) 



such that 



n* 



SK = - n ( MM - Mrus + mgrhs + nisarhsa + "^/i^^/i + "T-jjm-j} ) + njQ f m^m^ + msamf^ 
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(S53) 
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(S54) 



One can see that all the terms in Eq. ( 95^ are factorized with respect to the replica index 
a except the terms from the 3rd to the 5th line. In addition, we also have to linearize 
(iha)"^ and (ifja)'^, and decouple {iha){i'i)a), which will finally become the delta functions 
of ha and rja by integrating the corresponding ha and fj^ respectively. To achieve this we 



re-write n$ as 

n<l> = n{css-qss) + n{cssa<T-qs 



+iogm 
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(S55) 



where 



U 



and 



V 



' ^qss qgh ^^V qssa ' 

^sfj qfifj '^qfjfj qfjscr 

\qssa ^hscr qfjsa ^qssaa / 
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(S56) 



Va 






(S57) 



We can then adopt multivariate Gaussian integrals to linearize w'^ -lA • w and y^ ■ V ■ -jja, 
such that n$ becomes 



n$ = n{Css-qss) + n{Cssaa-qssaa) + log 



47rV|W| 
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(S58) 



We continue the calculation by: (i) expanding the exponential function in the 2nd line as a 
power of n, such that in the limit of small n one can make use of log[l + nC + 0(n^) + ...] ~ 
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nC to simplify the expression; (ii) collecting terms with factors of ih and ir), and the 
integration of h and 7] gives rise to the delta functions 



5[V+ (^3 + X2 + {Csf,-qsf,)s + {Cfisa " qfisa)s(T - m^j)] . 



Integrating h and 77 lead to 

1 



n 



X log 



X 



27rV|V| 
l-o- 1 + 0- 



dxidx2e 2^ ^ '^ 



4vrV|W| 

EE 



dzidz2dz^dzie '^' 



S = ±lcT=±l 



-e 






X exp 2;is + z^sa + m^s + nisaScr + (c^s^-gssajo- + 7- 



Finally, we sum over s and a to give 
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(S60) 



(S61) 
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(S62) 
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By using the above expression of $ and Eq. (9531). we can write lim„_>.o \1/sk as 

lim ^sK = -(mm - Mrus + msTJis + nisarrisa + rn-^rhf^ + mf^rhA + Jom^m^ + msarrifj 
+ (QQ - Qqss + QssQss + q.hQ.h + ^Is^sv + QssaQssa + qhh%h + ^;,^9fe^ 



j2 , , 

- Y (^ssg/,/, + (g,/,) + 2g/^^gss<x + 2g^,^gsr) + qi^ijQssaa + (g,}sa)^ ) 
j2 , , 

+ {Cssaa - qssaa) + / Vzlog I 2 COsh[zi - Z4 + ms-msa]e~''^'""'~^'""''^ + e'^^'""'~^'""''^ 

(S63) 
where J T>z represents the multivariate Gaussian integration 

r= I dzidzidzsdzie'^"^^'^^, (S64) 

and the functions Q± are given by 

Q±{z2,zs,mf^,mfj,c[,c) = ^== / dxidx2e~^'' "^ '"^ 

27?^/ |V| J 

such that q and c represent vectors of the variables of q and c respectively. 

2.1 Saddle point equations 

Since Sgx oc e^"*sK^ Q^g ^g^j^ evaluate 

— In ZQ,yi = — lim = lim \E'ot^ (S66) 

by the method of steepest descent such that \l/gj^ corresponds to the extremum of \E'sk with 
respect to M, Q and all the 48 variables of m, rh, q, q, c and c. Indeed, all the conjugate 
variable m, q and c can be extremized and expressed in terms of the variables of m, q and 
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■n 
(S65) 



c. To facilitate the presentation of the saddle point equations, we will denote the function 
in the logarithm of Eq. ( 963|) by T{z, m, q, c), i.e. 

J'(i', m, q, c) = < 2cosh[2:i-Z4+m5-m^^]e"^'^''""^'"""^ j^ ^(css.-qss.)+i 

X e-^i-^*-'^=-™-fi_(z2, ^3, m/,, m^, q, c) + e^i+^^+'^^+'^-fi+(z2, ^3, rn^, m^, q, c) 

(S67) 

We first differentiate the above expression with respect to M, Q and C and obtain the 
expected relations: 

m, = M, (S68) 

qss = Q, (S69) 

Css = 1. (S70) 

We then differentiate Eq. (S l63|) with respect to m^, mso-, '"^^ and m^^ to obtain the following 
relations which are identical to Eqs. (171)-(IT0l) of the Ising model 
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(S72) 
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(S73) 




rhfi = 


= JqITIscj. 




(S74) 


1 respect to the individual variables of q and c gives 
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Css = C+ 2 Chh 


(S75) 


Ih = J^^sh 






Kh = J^^'sh 


(S76) 


qsf, = J\hsa^ 
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(S80) 


ksa = J^9sf„ 
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Qf/sa •J Qfjsay 
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(S83) 


J' 






J' 

Cssaa r\ ^flV 


(S84) 
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The remaining tasks are to differentiate Eq. ( 9631) with respect to individual variables 
of rh, q and c, which involve differentiating the complicated function $ in Eq. (S l62|) . We 
first differentiate m^ and rhscr which give us 



m, = / Vz 



T{z, m, q, c) 



2sinh[2;i-Z4+ms-msa]e (^''»-^''»-) _^ g(£.s<T-<?s.a)+7 
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— e 
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(S85) 



X 



T{z, m, q, c) 



-21— 24— ms— 'nis 



2sinh[zi-Z4 + ms-msa-]e (^''»--9»=-) _^ e('5.»<T^ss<T)+7 



'fi_(2;2,2;3,'^ft,^»?,q, c) + e 



zi+zi+rhs+rhs 



'Q+{z2,Z3,rhf^,rhfj,q,c) 
(S86) 



We remark that m^ = M is known in Eq. ( 95^ and one should instead extract the value 
of Trig from Eq. (S IS^ . The differentiation of Eq. (SIHS]) with respect to m^ and rhfj involves 
differentiating the step functions in Q± in Eq. ( 965|l and should be taken with extra care. 
We first differentiate Q± with respect to rhf^ 
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(S87) 



where we arrive at the last line by integrating xi in the delta function, such that Xi in the 
final expression is substituted by xf = —Z2 ^{Csh~Qsh)~'r{(^hsa ~ %sa)'^^h- ^^^ ^^'^ further 
simplify the above expression using the properties of multivariate Gaussian distribution 
and the definition of error function. If we denote the element in the i-th row and j-th of 
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V by Vij, the above expression becomes 
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where erf (x) is the standard error function. Similarly, the differentiation of Q± by rhf^ is 
given by 
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(S89) 



/ 



where x^ = — -23 =F(csr}— '?sj7)=F(cr}so- — ^r)so-)+'^r)- Finally, we can differentiate Eq. ( 963|) with 
respect to m^ and m^j to obtain an expression for mf^ and ?7i^ in terms of dQ^/drhf^ and 
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(S90) 



32;i+24+ms+mso- + 

drhf, 



(S91) 



We continue to differentiate with respect to the q variables. One can make use of the 
following lemma to simplify the calculations. 

Lemma Given a symmetric m x m matrix U, if we denote Uij to be the element of U in 
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the i-th row and j-th column, then 
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1 d 



TfTii-X 



duij ^p\ 2^p\du. 



z'U-'z 



']^jvzf{z){~{U-%,+ [zi{lA-%; + Z2{U-%, + --- + z„,{pi-') 

Jvzf{z)\ - {U-% + [zi{U-')u + z,{U-%, + ■■■ + z^{U-% 

zi{U-%j + Z2{U-%j + ■ ■ ■ + z^{U-^\ 



X 



where we have made use of the relation 



OUij OUij 



1= ] 



(S92) 
(S93) 



for the second term in the curly brackets. For f{z) to be a constant, one can show that 
-^— J Vz = for all i and j by the above lemma. 

To continue the calculation, we make use of the above lemma and denote 



Di 



VnogTiz,m,Ci,c)\-{U-%+ Zi{U ^)tl + Z2{U~^)i2 + Z^{U ^)^Z + Zi{U )i4 



/-l^ 



I = J 



I Vz\ogT{z,mA,c)\ - {U-% + [zi{U-^\i + Z2{U-\2 + zz{U-%^ + Zi{U-^),i 



Zi{U-')ji + Z2iU-')j2 + Z3iU-')j3 + Z^{U-%i 



(S94) 



The differentiations of Eq. 



631) with respect to Css, Qss, Cssaa and Qssaa lead to 
c =1 



1, 

1 - 2D44, 



c =1 



q. 



(S95) 
(S96) 
(S97) 
(S98) 
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where Qss = Q is known and one should extract qss from the right hand side of Eq. (S l96|l . 
The factor 2 in Eqs. (S l96|) and (S l98|) comes from the fact that 2qss and 2gsso-o- are the 
elements of the matrix U. The differentiations of Eq. ( ^63l) with respect to Cssa and gssa 
lead to 



Vz 



T{z, m, q, c) 



2 cosh[zi- Z4+ms-msa]e (^'"'-^»''-) ^ Qii^ss,y-qss,y)+i 



X 



-^1— 24— m.s— rris 



'fi_(2;2,2;3,m^,m^,q,c) + e 



zi+Z4,+ms+ms 



n+{z2,Z3,mf^,m^,q,c) 



Qssu — -i-^14 + Csscr- 



(S99) 
(SlOO) 



From the definition Eq. (Sl44l) of g"^'^ for a = /3, the size of self-sustained clusters is indeed 
given by r = (1 + Cssa)/2. One can also differentiate Eq. ( 963|) with respect to 7 to show 
this relation. We then go on to differentiate with respect to c^^, c^^^, c^^ and Cfisa by noting 
that 



9fi± dn± dn± 



9c,h dksa 



T 



drhi 



"h 



(SlOl) 
(S102) 



The differentiations of Eq. (S l63|) with respect to c^^, g^^, c^^^ and g/^^^ lead to 



Csk = I T^^ 



^(Csso— <JsstT)+7 



F{z, m, q, c) 



— e 



— Zl-24-ms— TTlsCT 






+ e 



zi+Z4+rhs+rhsa + 



'^hsa ^sh 



%sa - --^24 + Cf^s„ 

Similarly, the differentiations with respect to Csfj, Qsri, Cf^sa and g^so- lead to 



(S103) 
(S104) 
(S105) 
(S106) 



■^s-q 



Vz 



J={z, m, q, c) 



^(Csso— gssCT)+7 



3— 21— 24— ms— msCT 



9fi_ 



(9c, 



+ e 



21 +24+rhs +msCT '^ + 



s»? 



dc 



Qsfj — ~-Di3 + Csf, 



sf) _ 

(S107) 
(S108) 
(S109) 
(SllO) 
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Finally, we differentiate Eq. (S l63|) with respect to c^/^, qj^f^, Cfff,, qffi), c^ and g/^-, which 
are elements of the covariance matrix V of the Gaussian distribution in VL±. We thus make 
use of the above lemma again and denote 



D^^ 



1 



47rV|V| 



rfxidx2e-^^^^^"^^ <j -(V-i).. + zi(V-i).i + Z2{V-%2 



Z2 + XI ±(c^^-g^^)±(c^^^ - %J-rnA Z3 + X2±(c,^-g,^)±(c^,, - g^,,) 



-m^ 



? = J 



^y(ixidx2e-^^"^^"'^^{-(V-i)., + [zi(V-i)a + ^2(V-i)i2][^i(V-i),i + ^2(V-i),2]} 



27rV|V| 

xe[(z2+xi±(c^^-g;^)±(c^,^ 



'"^/ij \ZZ + X2^{Csf,-qs'i))^{Cf,sa -qfjsa)-mfj 



(Sill) 



The differentiations of Eq. ( 963|) with respect to c^^, g^^, c^^^, gj^^^, c^^ and g^- lead to 



Cm = / T^z 



T{z, m, q, c) 



3(Cs3CT-<?SSct)+7 



2e 



-^1— 24— ms— m-s 



■■ n^- _]_ 2e-^i+^4+^s+rnso- 7~)^+ 



%ft = -2^22 + Cy^h 



VV i-^'2 



(S112) 
(S113) 



J'(£, m, q, c) 



^(Csso— qss<T)+7 



2e- 



—z\~Z4,-ms-rhsa rP'- _L o„2i+-Z4+m.s+m.sCT n^+ 



d;;- + 2e 



D 



22 



Qflfi = -2D33 + C 



vv 



'^^ /^V(.-',m,q,c) 



3{CsSCT-l3sSO-)+7 



-zi-Zi-'tfis-irisa -ry^- _|_ „2i+24+ms+nis<T r)f^+ 



'12 



g/,^ = -£"23 + cf^f, 



(S114) 
(S115) 



(S116) 
(S117) 



We iterate all the above equations numerically to obtain the solution of all the unknown 
variables. 
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2.2 Final expression of In 2^1 



SK 



Finally, we use the relations derived in the previous subsection to express In .Zsk in terms 
of M, Q, rhg, qss and all the other variables of m and q, which lead to 

1 J2 

— InZsK = -Mrhs - Joms^mf, + {Q - l)qss + -^(^vv ~ 1m) 

+ Y ( {ish) + '^%r^1ssa + '^Qhsalsv + QmQssaa + (qijsaf 
t2 

+ jDzlog <2cosh[zi~Z4 + 'ms — Jomfi]e~ ^'^hfi~^hfj) -^ ^ y'^'hfr^'hf,i^ 



e-^^-"^-™=--^o™'^fl-(z2, ^3, M, m,,, q, c) + e^i+^4+'^^+^°™^0+(z2, ^3, M, m,,, q, c) 



(S118) 



where J Pi* represents the Gaussian integration 

with the covariance matrix U given by 

/2qss/J^ q,h Qhsa (lhfj\ 

II ^ j2 Qsh Q Issa Qs-q 

^hsa Qssa Qssaa Qf/sa 



(S119) 



The functions Vt^ are given by 

fi±(z2,Z3,^,"^sa,q,c) 



2tiJ\V\ 



i^Z2+Xi ±J'^{c^f^-q^f^)±J'^{Csfi - qsfj)-JoMj 



X [Zs+X2±J {Cf^^^-qf^^J±J {Cf,sa-qf,sa)-Joms 



where the covariance matrix V is given by 

t2 / -'- ^; (^ssa Hssa 



V = r 



c.« 



Q.ssa J- Q's 



(S120) 



(S121) 



In the limit of Jq — ?> oo, which reduces to the Ising model case, one can show that: (i) 
VL± — )■ Q{Mmsa) when M 7^ and nisa 7^ 0; (ii) with the ansatz rhs oc Jo, cosh(2;i — 2:4 + 
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rhs — Jo^jy) ~ cosh^ifis — Jorrif^) and e='=('^i+^4+ms+Jom^) ~ g±(ms+Jomrj)^ which imphes InZsK 
is independent of the matrix U and each c and its corresponding q variable will become 
equal; (iii) these lead to rhs = PJqM which implies (ii) is a self consistent assumption. 
These results show that in the limit of Jo — )■ oo, Zsk in Eq. (HUE]) reduces to Rising in 
Eq. (9331) of the Ising model. 



2.3 Derivation of In 5(M, Q) 

We will derive an expression for the spin entropic contribution \n.B{{Ma},{Qai3}) given 
by 



lnS({M4,{Q„4) = lim- 

n— >-o n 



,s n 4^ - *'") n * (^^^ - <3"^ 



N 



(S122) 



We start with 












= Tr 


J 271 


a/3 


y 27r ' 


a 


vr lly 27r 

a/3 







(S123) 



where we have used the change of variables Mq, — )■ iNM^ and Qa/3 — ^ iNQap to arrive at 
the last line, with n\l/B given by 



n^B = - XI ^°^" - Yl Q»^Q^f^ + log ^ Tr exp 

a a,^ k 



y^ MaSo + X Qa/3SaS/3 



. (S124) 



Using the replica symmetric ansatz Mq, = M and M^ = M for all a, Qa/^ = Q and Qajs = Q 
for all a 7^ /3, and (5a/3 = 1 and Qap = C ioi al\ a = (3, we arrive at 

n^B = -nMM - n{n - l)QQ -nC + log { Tr exp 

{Sa} 



M^^s^ + gf^s^j +^((7-^) 

a \ a / 



1 



-nMM - n{n - 1)QQ - nQ + log <^ —^= / dze'^ JJ 



.s=±l 



'2QSa+MSa 



In the limit of n — )■ 0, \1/b is given by 

1 /■ .2 

lim ^B = -MM + QQ -Q + -^= I dze~- loj; 
"^0 V2^ 



2 cosh \j2Qz + M 



. (S125) 



(S126) 
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By the method of steepest descent, we differentiate Eq. ( 91261) with respect to M and Q 
to obtain 



M 



Q 



f^TX 



dze'^ tanh y2Qz + M 



dze~^ tanh y2Qz + M 



(S127) 



^2n 



dze~- tanh^ y2Qz + M 



(S128) 



which are identical to the equation of states in the original SK model when the unknown 
M = (3JoM and Q = /3^J^Q/2. Substituting this solution into Eq. ( 9126p . we obtain an 
expression for In B{M,Q), given by 



In 5 (M,Q) = lim^B 

n->0 



2 72 



(3' J 



■ -^JoM' + ^^^—Q{Q-l) + 



dze 2 log 



2cosh[(3J^yQz + /3JoM 

(S129) 
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